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ABSTRACT 


Implantable pulsed— Doppler ultrasonic flowmeter development has 
resulted in designs for application to the aortas of dogs and humans, 
and to human renal and coronary arteric... A figure of merit has been 
derived for each design, indicating the degree of its precision. 

An H- array design for transcutaneous observation of blood flow 
has been developed and tested Jin vitro . Two other simplified designs 
for the same purpose obviate the need to determine vessel orientation 
One of these will be devcl- od in the next time period. 

Techniques for intraoperative use and for implantation have had 
mixed success. While satisfactory on large vessels, higher ultra- 
sonic frequencies and alteration of transducer design are required 
for satisfactory operation of pulsed Doppler flowmeters with small 
vessels . 
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I . INTRODUCTION 

This semi-annual report covers the period from September 1, 1973 r 
through February 28, 1974. It discusses recent developments in work 
toward achievement of satisfactory blood flowmeters suitable for 
measuring flow in the renal rrtery, both with implantable devices and 
with external devices operating transcutaneous ly . The use of custom 
integrated circuits and ultrasound are essential to accomplish this 
task in the approach described here. 

Implantable devices are essential to permit unhampered observa- 
tion of blood flow to specific organs in awake, unanesthetized animals. 
Transcutaneous observations avoid the necessity for invasion of the 
body in order to measure the flow. Instruments available previously 
for flow measurement have been less accurate and reliable than those 
described here. 

Determination of renal blood flow in man is considered an important 
factor in understanding the mechanisms operative in producing the 
diuresis observed subsequent to weight lessne- s and bedrest simulation 
of weightlessness and, therefore, are important in understanding the 
effects of space flight on man. 

Analysis of several system tradeoff considers : ions has been com- 
pleted in the development of a pulsed- Doppler ultrasonic blood flow- 
meter. Included are pulse burst length, repetition frequency, filtering 
and ultrasonic frequency choice and transducer diameter. These are 
drawn together into an optimal design procedure and a study made of the 
performance limitations. 

To complement the design development work, implementation techniques 
are also being developed using dogs to test the effectiveness of the 
designs as a prelude to their application to people. 

II. TECHNICAL RESULTS 

A . Implantable Pulsed Flowmeter (R. W. Gill and J. D. Meindl) 

Work on the realization of the pulsed flowmeter implantable 
electronics in intograted-circu.it form is continuing. The final circuit 
configurations have been chosen on the basis of minimum power drain and 
reliable operation over a wide range of supply voltages. These circuits 

have now been laid out as integrated circuits, and the IC masks made. 

Emphasis is, at present, being placed on development of a suitable IC 
fabrication process, and it is expected that working integrated circuits 

will be fabricated in the near future. Work is also proceeding on the 
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wideband telemetry link required for tile pulsed flowmeter. This has 
proved to be a more difficult task than was originally envisioned. 

The following discussion is a reproduction of a relevant paper on 
optimal system design of the pulsed Doppler ultrasonic blood flowmeter 
presented by the authors at the IEEE 1973 Ultrasonics Symposium held at 
Monterey, California on November 5, 1973. The work presented here was 
also supported by the National Institutes of Health under PHS Research 
Grant No. 1 P01 GM 17940- 04. 
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OPTIMAL SYSTEM DESIGN OK THE PULSED DOPPLER ULTRASONIC BLOOD FJ/WMETtiR 

1 1 . W , (Jill and J.1). Mcindl 
SLa.ti.ford K1 octrunJ .■ ’ Laboratories 
Stanford, C'.A 94305 

A BSTRACT . The pulsed Doppler blood flowmeter Has demonstrated Its utility in a number 
of applications, both clinical and research. The future is certain to see increasing use 
of this type of flowmeter for several reasons, including its convenience and stability and 
the unique data it can produce. 

In the design of the flowmeter a number of tradeoffs must be made. This paper presents 
a systematic approach to these tradeoffs, leading to a design procedure which optimizes flow- 
meter perlon .nee for a given set of physiological constraints. Sample designs for typical 
applications are presented. The particular flowmeter configuration considered is the single- 
transducer type with the transducer adjacent to the vessel. 


This paper deals with the tradeoffs involved in 
the design of the pulsed Doppler ultrasonic blood 
flowmeter, and the .selection of the various parameters 
of the flowmeter to optimize its performance subject 
to the physiologic^ constraints of the application. 
Specifically, the analysis is applicable to the case 
of a no»-di reel ional^ or quadrature directional 3 
pulsed flowmeter, with a single transmit/receive trans- 
ducer adjacent to the vessel. 'Hie analysis cm be 
modified to deal with an offset directional flowmeter 3 , 
transcutaneous measurement with or without a lens, and 
the case here separate transmitting and receiving 
transduc. s are used. 

Optimization of the flowmeter performance is taken 
to mean obtaining maximum range resolution and flow 
measurement accuracy suijject to the needs for unambig- 
uous measurement of range and velocity. The maximiza- 
tion of the Doppler signal-to-noise ratio for a given 
transmitted acoustic power is also an important con- 
sideration. * t 

A simple first-order theory is developed, since 
it is desirable to study the design tradeoffs analy- 
tically. Thus, for example, the beam pattern c>f the 
transducer is modelled by r theory which may deviate 
significantly fro..i the actual pattern in any given 
case, and the samp: volume of the flowmeter has to 

be somewhat arbitrarily defined, since the intensity 
falls off slowly at its edgi s. In addition, the 
problems of optimum demodulation of the Doppler signal 
have been omitted; this subject deserves considerable 
attention in its own right. 

The block diagram of an implantable pulsed flow- 
meter is shown in Fig. 1, and some typical waveforms 
in Fig. 2. 4 Only one demodulator channel is shown for 
convenience. A non-implantable version of this is 
obtained by eliminating the telemetry link and connec- 
ting the video and sync signals directly to the demod- 
ulator. Although other pulsed flowmeters differ from 
this in several respects, the following analysis should 
still apply to them directly. The parameters which 
must be determined are; ultrasonic frequency, burst 
repetition frequency, transducer diameter, transducer 
angle, burst length, rf amplifier and video amplifier 
bandwidth, arid audio filter cutoff frequencies. The 
factors influencing the choice of each parameter are 
discussed, and then a design procedure presented which 
combines all these considerations to give an optimized 
flowmeter design. The flowmeter performance and the 
effects of physiological and transducer constraints 
on it are then discussed. Finally, some typical 
designs are presented , 
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Fig. 1. Pulsed Flowmeter Block Diagram 



Fig. 2. Typical F) ov. meter Waveforms 
burst Length 

Fig. 3 shows tire beam and vessel geometry and the 
definition of some terms. The flowmeter resolution 
(X) is given by the projection of the sample volume 
(the volume of blood from which significant scattering 
is received at any one Instant) on the vessel diameter. 
Jt can be expressed as 

X =1/2 c T sin 0 + w COS 0 (1) 

where c-1.57 x 10 cm/sec is the velocity of sound in 
blood, T off is the effective burst length and w is the 
beam width. To maximize resolution T^.^- must clearly 
be minimized. It is determined by the' length of the 
drive burst applied to the transducer, ui.u by the band- 
widths oj the transducer am! nil the following elec- 
tronics tip to the sampl o/boid gates. In practice 









i'clf can be* defined ;>b the interval between the ~10dB 
points of Uie return i roiii a single reflector, measured 
at tlie output of the video ainpli i 1 er . 


be measured, according to the Doppler equation 

X , = 2 - \\ cmQi 
d c 


( 5 ) 
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Fig. 3. Beam and Vessel Geometry 
The bandwidth of the transducer is given ' : f /q, 

- ei 

where is the transducer Q, and 1% is the ultrasonic 
frequency, The rf amplifier bandwidth is set equal to 
the transducer b.v Iwidth; a narrower rf bandwidth would 
degrade resolution, while a broader one would unnecessar- 
ily degrade the si; igil-to-noise ratio. 5 Since the mixer 
folds the spectrum in half about zero frequency, the 
video bandwidth is one half the rf bandwidth. The 
length of the drive burst is made equal tn the recipro- 
cal oi the transducer bandwidth 

T = Q t /f e (2) 

since a shorter burst would reduce the power received 
without significantly improving resolution while a 
longer one would degrade the resolution. The effective 
burst, length can then be written as 


where is the Doppler shift, v is the blood velocity, 
c is the velocity of sound anil f„ is the ultrasonic 
frequency . 

Burst Repot i lion Frequency 

The above discuss; ion implies that a high repeti- 
tion fre :ency is desired to allow the measurement of 
high velocities. However, the need to measure range 
unambiguously limits the repetition period to the time 
taken by sound to t r; vel to the opposite wall of the 
vessel and back to the transducer. Thus the repetition 
rate chosen is the highest one compatible with this 
latter constraint: 


_ c s'in 0 

1* ~ 2 D 

max 


( 6 ) 


where D max is the maximum vessel di; meter expected. 
Then, from ("1), the maximum allowable Doppler frequency 
is 


dmax 


0.225 


c Ai 9 
D 

max 


(7) 


Transducer Angle 


consider a jj.cwmci.ei 






velocity and volume flow. The velocity 
according to the Doppler equation (5) 
c f 

v = 


>asure both 
is calculated 


( 8 ) 


To calculate volume flow, the velocity distribution 
across the vessel in spatially averaged and multiplied 
by the vessel area, determined from the range of the 
near and far walls. Thus 


— 1_ 

v casd 

and D = l/2 c t 0 ~ Sil 0 


(9) 

( 10 ) 


f eff Q eff^ fu (3) 

where () ^ is the transducer Q modified to account for 
the effect of the rf and video amplifier bandwidth?.. 

Thus resolution is best for a low Q transducer 
and high ultrasonic frequency. Typically Q-j- may be 
reduced to about 5 by using a suitable piezoelectric 
matt : ial and acoustic loading and matching layers, 
giving Q Qn ~ 7. 

Audio Filter 

Kcflcctions from vessel walls and other body 
structures produce large low-frequency Doppler signals 
which tend to obscure the blood flow information. For 
this reason a lower cutoff frequency of 100 or 200 Hz 
is typically chosen. The* upper cutoff is chosen to be 
Just above the highest Doppler signal expected, and 
just below one half of the burst repetition rate. This 
is because the sampling theorem of comint in cat ions theory 
proves that the highest Doppler frequency which can bo 
nnambiguously measured is one half of the sampling rate;** 
in practice, filtering requirements reduce the maximum 
permissible Doppler frequency to about 90% of this value: 


where v is the spatially averaged velocity, t is the 
time between near and far wall echos, and D is the 
vessel diameter. The volume flow is then given by 

q = 1/4* '■' ! V~jjjjg 6 OD 

Differentiating (8) and (11) with respect to 0 and 
normalizing, the sensitivity of velocity and volume 
flow measurements to angle errors are given respectively 

by 


1 

dv 

= tan 0 

(12) 

V 

d0 



1 

dq 

= tan Q t 2 cot 0 

(13) 

q 

d0 




These are plotted in Fig. 4. For velocity measurement 
0 should be small for maxi mufti accuracy, but for volume 
flow measurement the sensitivity to angle error is 
lowest near 00°. Since angle or, irs of 1° or 2° seem 
unavoidable, it is essential that 0 bo chosen to 
minimize sensitivity Lo such errors. 


f. - 0.45 f (-1) 

dmax r 


where f|. is Uie burst repetition rate. Thus, for a 
given ultrasonic frequency and t ransdueer angle, a 
limit is set by i' r on the maximum velocity which can 
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Fig. 4. Velocity and Flow Sensitive, os 
to Angle Error 



Other factors which should also be considered in 
selecting 6 are-. sample size, maximum range fo; a 
given vessel diameter, and the distance traveller! aloir 
the vessel by the beam as it travels across the vessel. 
Th.i last factor determines the amount of relatively 
st) iglit vessel needed to make accurate measurements, 
wh j' i should be minimized, implying that a value of 0 
nea" 90 ° is desirable. The sample size is not strongly 
dependent on 9 for typical burst lengths and transducer 
size, but it is minimized for Q - 90 °. For a given 
vessel diameter, the maximum range limits the repetition 
frequency and so the Doppler frequency which can be 
measured, as given by (7). Again, the range is minimized 
for 0 = 90 v , and so a v« • c of 0 near 90 ° is desirable. 
The choice of 0 is a compromise between all these 
factors, with the need to minimize the sensitivity to 
angle errors the dominant consideration. For a flow- 
meter which is used to i. sure both velocity and volume 
flow, a value of 60° is i bus normally chosen. 


Ultrasonic Frequen cy 

As indicated earlier, and as would be expected 
intuitively, the ultrasonic frequency should be as high 
as possible to optimize resolution. However, for a 
given vessel size there is an upper limit on the Doppler 
frequency given by (7); for a given velocity this sets 
an tipper limit on the ultrasonic frequency, according 
to (5) . Equating the two upper limits 

sin G 


drnax 


= 2 - 
2 


f 0 tvrO - 0.225 


D 


i . e . f „ - - 


c tan 0 


8 . 9v 1.) 

ma x max 


(14) 


Using this value of the ultrasonic .frequency will then 
maximize the resolution subject to the repetition rate 
limitation on the Doppler frequency. Equation (14) is 
plotted in Fig. 5 for a number of vessel diameters and 
maximum velocities. 


Fig. 5. Frequency for Optimum Resolution 
as a Function of Vessel Diameter 
and Maximum Rlood Velocity 


One further consideration which must be discussed 
is the variation of signal-to-noise (s/n) ratio with 
ultrasonic frequency. It will be shown that this is 
not a critical factor in the choice of f 0 . 

Consider the scattering from blood at range r. 
Assume that all the t r ..remitted power is uniformly 
distributed in a beam of width w, area (l/4)nw 4 '. Then 


the incident power per unit area is given by 

\ ft 

a/4)rt~"w 2 e 


p. = 


(15) 


where P is the electrical power delivered to the trans- 
ducer, ij is the transducer transmitting efficiency, and 
a is the 1 attenuation constant of the intervening medium. 
For a single target with scattering cross-section A, 
the scattered power is just A.F.. If the scattering 
proco is first order and independent, as experimental 
results indicate, 7 the total scattered power is the sum 
of the powers duo to the individual scatters. Thus, 
it is given by A.P, multiplied by' the number of 
scatterers in the sampi volume. Ignoring the spreading 

of the burst in time ' ' .a. suming T - T) , the 

eli 

scattered power is th 

P s = [<l/4 n w 2 ) (l/2 c T) A.T^ 

= l/2 c Tn 0 A T) P e“ a ' r (16) 

vising (15), where n 0 is the number of scattercrs per 
unit volume. Since the scattering particles are much 
smaller than tile ultrasonic wave'ength ( r ' j 7|j compared 
to ^ 0.3mm) the scattering is governed by Rayleigh 
statistics, und so it may be assumed to be isotropic, 
with a scattering cross-section given by 


A 


(17) 


where A is a constant. 

The fraction of the scattered power which is 
received by the transducer is then given by the ratio 
of the transducer area to the area of a sphere of 
radius r. 
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where* is the transducer receiving el'Cio 1 ern-y , com- 
bin ing equations (l.O r (17) and (IK) gives 

2 •! -2a r 

P - 1) c. T n„ A,, tl 1> J„ c • (19) 

r ~r t 


where 7] - T| T)^ is tJie overall transducer efficiency. 

To make tho dependence of the signal -to-noi sc 
ratio on f a explicit, the following relationships are 
used : 

(i) .assume Rayleigh scattering, equation (17) 

(ii) the attenuation constant (a) is known to 
be approximately proportions, to frequency 
in the frequency range of interest 9 , so 



a -So 


f 0 


( 20 ) 


where a c is a constant , 

(iii) the burst length is given by (2) 

(iv) the receiver noise referred to the input can 
be assumed to be proportional to the receiver 
bandwidth, i.e. 


N 

r 


N„ ip 

% 


( 21 ) 


where N n is the noise per unit bandwidth . 

The s/N ratio of the Doppler signal is found by- 
substituting (2) and (20) into (19) and dividing the 
resulting received power by (21). 



Tl c Q “ n 0 A n d 2 -2a r 0 r 

L P t e (22) 

22 No r 2 


The electrical power (P^.) is usually fixed by 
electro. -iic and/or safety requirements, and the trans- 
ducer size is determined by resolution requirements. 

If we assume the transducer Q and efficiency are inde- 
pendent of frequency, then v.e have 



r 


-2a t i' n r 

c 


(23) 


Differentiating this with respect to frequency and 
setting the derivative equal to zero, the frequency 
which maximizes the S/N ratio is found to be 



This is plotted in Fig. 6 where r is assumed to 
be the range to the 'ar wall of the vessel, and a D is 
the weighted mean of the attenuation constant for the 
vessel walls (4dI3/cm/MHz) and 1 < uod (0 . 25dB/cm/MHz) 
assuming a ratio of wall thickness to vessel diameter 
of 5%. 

The shatlec area in Fig. 6 indicates the range of 
frequencies for a given distance (r) for which the 
.s/n ratio is within 3d]J of its maximum value. Tin s 
indicates that a range of frequencies of at least 4:1 
Is permissible without seriously affecting tho s/<N 
ratio. Thus the maximization of s/h ratio is a secon- 
dary consideration in the choice of ultrasonic frequency. 


Fig. G. Optimum Frequency for 
Maximum s/N ratio 

T ransducer Diameter 

Fig. 7(a) shows the beam envelope for a circular 
transmitting transducer of diameter d, assuming that 
it acts as an ideal piston.® The axial intensity dis- 
tribution is shown in Fig. 7(b). Notice that there 
are a series of maxima and zeros, their exact number 
being determined by the ratio d/X where X is the wave- 
length; d/A =10 for Fig. 7. 7'he point at which the 
last axial maxima:.; occurs is conventionally taken as 
the transition poi -t between the parallel-beam "near 
field" and the uni >rmly diverging "far field" where 
the axial intensity falls off inversely with range. 

This last maximum occurs approximately at the point 

r =0.25 d 2 /). (25) 

and this has been used as the definition of the maxi- 
mum usable range of the transducer. 

however, a number of comments are in order. First, 
the definition of the beam discussed above defines the 
boundaries as those points at which the intensity falls 
to zero. A more realistic definition would seem to 
be the* locus of points at which the beam po • or falls to 
a certain fraction of its maximum, since the flowmeter 
demodulator normally has some threshold power level 
below which it does not detect signal. ^ For example. 
Fig. 7(c) she.-, the ~GdB beam envelope, again usi- g the 
"piston" theory. 1 -^ Note that when the receiving sensi- 
tivity pattern is also taken into account, the returns 
from these points are actually 12dlJ down. The beam 
width at range rj_ is now less than l/2 d , and the angle 
of divergence of tho beam in the far field is about 
60% of that for the beam of Fig. 7(a). This definition 
of the beam is not only more optimistic, but also more 
const? ont with the observations of experimenters in 
the field of ultrasonic imaging.^ -2 

In addition to the above considerations, bokker 
et. al . have demonstrated that the uniform piston 
assumption is not very realistic, and instead have 
investigated tho eases of a transducer with simply- 
supported and clamped edges . 1 d j n noth cases the 
transducer vibration falls to zero at t:hc edges, con- 
sistent with observations of actual transducers. Fig. 
7(d) shows the axial intensity distribution for the 
simply-supported transducer for d/\=10; notice that 
the intensity variation in the near field is much less 
than for Die piston theory, and that the last axial 
maximum occurs at: about i*j= 0,11) d 2 /X . Analysis of 
the far field indicates that the angle of divergence 
of the ~6dB en. elope is 


rp = O.ftH X/d 


(26) 



(o) 


rtf 

L« 


i 

'*i>( O 

AXIAL UAfcAALjM 


(*> 


it) 




( 5 ) Tho burst length and bandwidths are chosen as 
discussed above. 



Fig, 8. Flowmeter Resolution as a 
Function of Frequency and 
Transducer Diameter 


Fig. 7. Transducer Beam Patterns and 
Axial Intensity Distribution 


Defining the maximum usable range of the transducer 
as the distance at which the -6dB points are 1.5d apart 
(i.e. the beam width is degraded 50'i) , equation (2G) 
gives 

r =0.85 d />, 
max * 

= 0.85 d 2 f„/c (27) 


Thus, the transducer is quite useful at ranges up to 
three times the near field distance given b (25). For 
a given application, the resolution is maximized if this 
maximum usable rang; is set equal to the distance to the 
far vessel wall; i.e. 


r 

max 

i.e. 


D 2 

max 0.85 < 1 f ,, 

sin 0 ~ c 


d 


■ D -j 1/2 

1.18 c max ' 

. f a S3 e J 


(28) 


Optimal Design Procedure 


The above considerations can now be drawn together 
to provide a logical sequence of design stops which 
optimize the flowmeti r performance within the limita- 
tions of the physiological situation and transducer 
parameters. The steps are described below: 


(1) 

( 2 ) 

(3) 


(4) 


Choose G, the transducer angle. Normally 
G =60*. 


From the maximum vessel -size expected, the 
repetition- rate is calculated from (6), 
ensuring unambiguous range measurement. 

'Fliis sets a limit on the highest Doppler 
signal, according to uq. (7). Choose the 
hi; best ultrasonic frequency compatible with 
this, according to eq . (Id) or Fig. 5. This 
ensures maximum rosolul ion subject to the 
need for unambiguous > Jocity measurement. 
Check to see if tho S/ ;i rat i <> is adequate. 
Choose the smallest transducer compatible 
with the maximum range required. Sir- ituting 
cq, (14) into (28) given 


r- 1 — 

(10 ii v 
I max max 

‘ |c .till 0 tan 0 


1/3 


( 29 ) 


Performance Limitations 

The flowmeter resolution is plotted in Fig. 8 as 
a function of frequency and transducer size, assuming 
a parallel beam of width d for simplicity. Notice, 
however, that the design procedure allows the ultra- 
sonic frequency to be expry. ed in terms of v max and 
D max , the physiologic, 1 con: "aints, according to eq . 

14. When this is taken into ecount, the resolution 
can be written as 

X = D P.U Q ,. -S££«ye) + cot 0 (10 — ~ COS 6) 1 ^ 
max r eff c c ^ 

(30) 

Thus !'■ : ratio D max /X, roc. hly speaking the number of 
indepei dent samples icsol •. . ..ble across the vessel, is a 
function only of the maxi- urn velocity expected, the 
transducer Q, and the tra; wducei* angle. This function 
is plotted in Fig. 9 for G_G0° . Together with tho 
an; 1 e-sensitivity curves of Fig. 4, this is probably 
tin most important fun unental limit on flowmeter per- 
form:!' ce. Other limits are either built into the design 
(the , -.ximum velocity and ve;: 1 diameter) or not fun- 

damental (the s/n ratio can be improved by increasing 
the transmitted power or improving the processing cir- 



Fig. 9. Number of Samples Across Vessel as a Function 
of Maximum Blood Velocity. 
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Table I presents some sample T] owmeter design:- J’ur 
a variety of applications. It is interesting to note 
that the designs for the large vessels arc very .similar 
to those alia - 4y used, and also that (he figure oi merit 
(the number oi independent samples) is quite reasonable. 

In summary, the ('actors Influencing the choice oi' 
the major flowmeter parameters have been discussed, and 
a design procedure developed which optimizes the flow- 
meter performance for a given physiological situation. 

The flowmeter performance is then specified, and its 
dependence on the physiological, parameters and the 
transducer properties is explicitly stated. 
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Ta bic I 

sample design s 


VESSEL 

dog 

aorta 

huma n 
aorta 

MAX. 

DIAMETER 

2.0 

4.2 

MAX . 

VELOCITY 

200 

150 

TRANSDUCER 

ANGLE 

60 

60 

BURST REP. 
FREQ. 

34 

16,2 

ULTRASONIC 

FREQ. 

12.1 

7.7 

TRANSDUCER 

D1AM. 

1 .8 

3.4 

SAMPLE SIZE 

1.7 

2.9 

FIGURE OF 
MERIT 

11.7 

14.4 

OPTIMUM f 
FOR S/N C 

4.3 

2.2 


liu man human Units 

re n . art . cor. art, 

0.9 0.3 cm 

65 200 cm/ sec 

GO 60 degrees 

76 226 kHz. 

83 81 MHz 

0.5 0.3 ram 

0.4 0.3 ram 

25.6 11.7 

9.6 29 MHz 
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B . Transcn taneous FI c wmotcrs 

1. Arterial Visualisati on System (C.F. Hottinger) 

The goal of the transcutaneous blood flowmeter project is 
development of an ultrasonic system capable of quantitative measurement 
of instantaneous flow through the renal artery. To perform this 
measurement, the transcutaneous system must be capable of real-time 
visualization of the effective vessel lumen. The arterial visualization 
syst- in cur. ently under development generates one or more Doppler B-scans. 
Sensitive only to moving targets, each scan displays the projection of 
the effective lumen area onto a reference plane. 

An appropriate choice of reference planes permits an unambiguous 
determination of eft* olive vessel size and orientation. For this 
reason an H- shaped conf igura tion. of linear transducer arrays has been 
developed; when the H- array is aligned over the vessel axis, the vessel 
size and orientation can b*. estimated reliably. A separa i transducer 
element, not part of the visualization system, can then be used wit’ 
any of several pulsed Doppler demodulator schemes to make a reliable 
quantitative measurement of flow in real time. 

The H-array configuration has been implemented in an assembly 
involving 78 transducer elements, with center- to- center element spacing 
as close as 1 mm. Element s; rcings in this range give the type of 
near- field resolution acceptable with the lai’ger arteries; closer 
spacing will achieve the resolution appropriate to the renal artery. 

In addition, high voltage multiplex circuitry was introduced to permit 
the higher transmit levels needed for transcutaneous scanning. 

The Doppler B-scans indicate the points under the arrays returning 
total Doppler- shi f ted audio power above a pre-set threshold. Even 
when this f low- veloci ty vector has no component in the direction of 
the scanning transducers, Doppler shifted audio power is detected. 

Such mechanisms as beam spread, finite transit time, and flow turbulence 
contribute to the apparc- t Doppler shift. 

The low power levels involved, however, prompted consideration of 
other array configurations that avoid perpendicular scanning of the 
velocity vector. Two new configurations, differing radically from the 
H-array, have been developed. Both new systems combine a first- moment 
processor (discussed below) with the Doppler imaging processor in a 
new configuration to simplify greatly and Increase 1:he accuracy of 



10 


ultrasonic volume How measurement. Described as performing a 
"canonical" measurement of the flow, the method consists of breaking 
up an arbitrary plane of incremental areas. The incremental flows 
across these areas are summed to find the entire volume flow. There- 
fore, no assumptions re;; rding orientation, turbulence, or velocity 
profile are necessary to estimate reliably tile total flow. 

The significance of the canonical method lies in its solution, in 
concept, of the problem of flow measurement for a large number of the 
vessels of interest in the body, in the coming year, the simpler of 
the two canonical methods will be implemented. involving a single 
linear array for visualization, as well as a uniform illumination 
transducer for first-moment detection., the technique promises to be 
useful both for transcutaneous and intraoperative measurements . 

Initial tests with arterial models and in vitro specimens will provide 
a basis for meaningful in vivo measurement i later. 

Further details of this development are given in the following 
reproduction of a pap- • presented at the 1973 IEEE Ultrasonics Sym- 
posium held at Monterey, California on November 5, 1973. 



11 


AN U LULA 5 ON IC SCANNING SYSTEM FOR ARTERIAL IMAGING 

C.F. Hot i;ing-=r nd J.D, Kcindl 
Department of Electrical Engineering 
Stanford University 
Stanford, California, 94305 

ABSTRACT. A transcutaneous system fo>* real-time imaging of arteries is described. A 
transducer probe containing three linear arrays of elements scans tissues within the element 
Fresnel zones. Each element is multiplexed in turn to a pulsed Doppler processor to permit 
real-time visualizations of movie . targets. Applications of the motion-indicating B-scans 
are disc -ssed, with special attention to problems in transcutaneous measurement of blood flow. 


A transcutaneous Doppler imaging system has been 
developed to disolay real-time visualizations of some 
of the larger peripheral arteries and veins. The 
probe for this ultrasonic imaging system consists of 
an array of pulsed transducer elements, each scanning 
tissues within be elc-r-nt Fresnel zones. A pulsed 
Doppler process-- * is demultiplexed to each of the 
transducer clem ts in turn to achieve real-time 
displays of mo v g 1. -gets and suppression of sta- 
tionary target di.'-j .*= . 

A convenient j al-time technique for visualization 
of the larger arteries and veins by transcutaneous 
moans would be valuable both as a first step in trans- 
cutaneous blood flow measurements and as a procedure 
by itself performed to evaluate the condition of 
vessels. Accurate quantitative blood flow measure- 
ments could be performed with a suitable transcuta- 
neous ultrasonic flowmeter since vessel location, size, 
and orientation could bo reliably determined with the 
visualization system, by itself, the imaging tech- 
nique would aid in the diagnosis o.f several vascular 
disorders since lumen size could easily be measured to 
iOCiu,s constrict ions uxuiig <,uo vessel length , In 
addH 'on, a system with adequate resolution would 
allow changes In the diameters of arteries during tb-? 
cardiac cycle to be wear red. From these changes an 
estimate of vessel disl sibility could be derived. 

Transcutaneous visualizations of vessels using 
ultrasonic Doppler tt-chniqo , have been reported in the 
literature, but these opera' .as involved multiple 
scans by a single, transduce r transducer pair. Other 
work'* 1 with a sirc.,le pulsc-c* ■ transducer generated 
real-time A-scans across vessels, but relied on specu- 
lar reflections off the vessel wa l 1 s to determine 
vessel size. Not only are such specular reflections 
inconveniently sensitive to transducer alignment, but 
the blood-tissue interface within diseased vessels may 
be difficult to locate using simple pulse-echo tech- 
niques under even optimal conditions. Methods utiliz- 
ing Doppler techniques rosj <nd only to moving targets, 
so delineation of the volumes actually containing flow 
is possible. However, these scans should be accom- 
plished in real time to permit convenient examination 
and reduce susceptibility to motion artifact. 

The physical constraints that allow the arterial 
scanner to ; lay a crucial role in the transcutaneous 
measurement of blood flow are determined by suL titut- 
ing the Doppler relation into the equation for volume 
flow, Q (in ml/soc) through a lumen of area A. 

5 ■ ( tt ) t5Te <^> A <» 

where c/f 0 the .ultrasonic wavelength, 6 is the 
angle between the average flow vector and the ultra- 
sonic beam path, and (Af) iu the spatial average 
Doppler sli li t measured across the lmwn. Of t he four 
factors involved in Q, only the first, ultrasonic 


wavelength, is known a priori. Tbo average Doppler 
shift (Af) Klusi be measured quant j tati vely using an 
appropriate Doppler flowmeter, while the direction 
cosine and area can be determined by a visualization 
of the vessel volume containing flow. 

Figure 1 shows the transducer array configuration 
that was developed to use the arterial scanning system 
in the transcutaneou_ mcasu ement of blood flow. 


ARRAY £ 2 



Figure 1. H-array of Transducers . 

Three linear arrays of transduce.* elements are combined 
to form an U-shaped probe that is placed flat against 
the skin of the vessel that is b.ing examined. The 
probe is me.' at to be aligned so th t the cross-bar of 
the H is directly over the projection of the vessel 
axis on the skin surface. Each of the transducer 
elements that are shown in the figure to be flat 
against the skin are part of th- arterial imaging 
system. The range-gated processor examines a number 
of depth increments under each transducer element in 
order to make a qualitative .judgement as to whether 
flow exists at these points or not. The single trans- 
ducer element that is shown in the figure to be tilted 
away from the skin surface is not part of the imaging 
system per se but is used to make a quantitative 
measurement of the average Doppler shift across the 
lumen . 

Figure 2 illustrates the orthogonal real-time 
displays of the vessel that are generated by the* arte- 
rial scanning system when the probe is properly aligned. 
The parallel segments of the H-array generate scans o [ 
the vessel cross section; when each image :s centered 
on the row of I is display, the pro bo is properly 

aligiu- , and Urn lumen area can be estimated from these 
cross-suet ional images. The orientation of the length- 
wise* scan image , generated by array S3, the cross-bar 
of the II, indicates the angle of the blood velocity 
vector with run peel i:< the skin. In this way the 
arterial scanner measures two of the three quant i t les 
that must bo measured to estimate volume blood flow, 
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Figure 2, Visualizations Generated by H-arruy. 


Figure 3 shows a simple block diagram of the 
processor associated with each one of the 11: linear 

arrays in the H~ shaped probe. Tire simplest j -ccxror 
utilizing this scheme uses linear audio filters and is 
not capable of distinguishing flow direction. A later 
version utilizes an off-set frequency technique and 
zero-crossing ck vector to allow determination of flow 
direction. 



Each transducer element is activated ior 10msec 
before the mult j plexor is slopped if Doppler shifts as 
low as 100 If?, are to be detected. The scope trace Is 
then a-lso stepped in Ibo y-cliroction every 10 msec; 
with 10 transducer elements in an array wo then achieve 
10 frames per second . 

The later processing system capable of distinguish- 
ing flow direction operates similarly. However, in 
this case the reflected signal is mixed with a fi - 
quency offset by 1 kHz from the drive frequency. Also, 
instead of a linear audio filter bank, a bank of zero- 
crossing detectors and counters is used. During the 
10 msec a particular transducer is being pulsed, 
counters indicating less than 30 counts indicate flow 
away from ho transducer, and counts greater than 'll 
Indicate flow toward the transducer. 

The first JI~ array to be assembled contains a total 
of 39 cle. ... nt pairs. Each clement pair is in turn 
multiplexed to a hybrid transformer scheme allowing 
common-mode drive and differential receive configura- 
tions. The element pair dimensions in the cross-arrays 
are 1.5 tan x 3 imn, while those in the lengthwise array 
are 2 ram x 3 mia. The transducer material is I.TZ-2, 
utilizing an energy ab orptive backing and quarter-wave 
transformer layer to increase transducer bandwidth. 

The transducer multiplexer now used is assembled 
from commercial JFET components capable of switching 
20 volts peak to peak; when these are replaced bv 
higher voltage devices, more power can be delivered to 
the far end of the transducer Fresnel zones. This 
capability will allow peno ration of the layers of fat 
and muscle that cover vessels like the carotid arterv 
system. in this way routine real-time scans of these 
and other important arteries can be made, and trans- 
cutaneous measurements can be made of the flow through 
these vessels. 
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Figure 3. block Diagram of Doppler Imaging Processor. 

The simpler system, using linear audio filters, 
operates as follows. Every 50 microseconds the driver 
generates a 1 gsec burst of 6 MHz sinusoid. This 
signal passes through Lbc* transducer multiplexer to a 
single one of the transducer elements. The r< fleeted 
ultrasonic signal, detected by the same element, 
returns through the transducer multiplexer to the 
Doppler receiver. Here the returning signal is mixed 
with the local or cilia tor output, and the resulting 
video is range gated through a demultiplexer to an 
audio filter bank. Each of these audio bandpass filters 
is used to detect u Doppler shift, d reflection from a 
particular depth under the transducer currently being 
pulsed , while the level of the shift detected modulates 
the CUT intonsi! v. At any given time the state of the 
audio filter mu Li i plexor determines the x-displacement 
of the CRT trace. Since the entire audio filter bank 
is sequenced beginning with each transmitted pulse, the 
scope trace covers the entire horizontal axis every 
50 psoc. Likewise at any given time the state of the 
transducer multiplexer <la tormin'-s the y-rii. a placement of 
the trace. 
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2, Doppler Spect rum Fi rst Mom . -» t Detect or (Levy Gerzberg) 

To measure blood flow, the average velocity of blood particles 
in the cross-section area, the lumen, of the vessel needs to be measured. 

It can be shown that, in the case of uniform illumination, of the blood 
vessel with the ultrasonic beam, the average velocity is proportional 
to the normalized first moment of the Doppler power spectrum of the 
returned signal. The - type first moment detector was found te be 
the optimal dt <odulator for use with the various types of Doppler flow- 
meters. This detector was first simulated on a digital computer for 
optimal system design and then was built with commercially available 
discrete components. With the aid of computer-aided design (CAD) many 
more alternatives could be considered titan would have been possible if 
each had had to be physically constructed. The accuracy of the otector 
was found to be: 2% with a single- frequency input signal and 5% with a 
wide- band signal of 18 dB s ignal- to- noise ratio. In the future a 
monolithic form of the first- moment detector will be developed. 

As a first step in the development of the pionolithic form of this 
detector, CAD has also been used. The computer is being used both for 
the rimulat on of models for some of the de\ ices used in the detector 
and for optimal design of the integrated circuit when processing 
limitations and parasitic effects ax^e taken into account. 
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C. Implementation To- ’'.liquej (W. R. Freund) 

Tlie study to d- termine a valid calibration scheme 1'or measurement of 
blood flow in the renal and iliac arteries has continued. Telemeli red 
information ha been collected from eleven animals during periods of 
rest and excercise. These data will be summarized and analyzed to 
establish the normal range and velocities which can be exp ited for renal 
and iliac flows in unanesthetized canine: . Arteriograms have been com- 
pleted on four of these animals to measure the vessel lumen area and the 
amount of occlusion produced b 3 ? the mall- vessel CW cuffs. The diameters 
of the three left and one right renal artery range in value from 2.5 to 
3.0 mm with a mean diameter of 2.75 mm. The diam Lers of the four iliac 
arteries ranged in value from 2.9 to 4.3 mm with a mean diameter of 3.55 mm. 
In one of the animals a significant occlusion was produced by an additional 
CW cuff placed on the right renal artery. The lumen was measured to be 
less than one mm in diameter as a result. While the kidney a -eared to 
be well perfused and free from infarcts, pressure cathether measure- 
ments showed t at the an!.- 1 was sevei fy hyper b nsive . A direct 
bleedout calibration produced maximum volume flows oC SO ml/min (approx- 
imately l/3 the expected flow for the 20 kg animal) . Post-mortem exam- 
ination revealed that the occlusion was due to an incorrect choice of 
cuff size rather to exc-.ss.ive formation of scar tissue. With the excep- 
tion of this case, no vessel occlusion has been observed. 

The development of special purpose measurement techniques has con- 
tinued. Three specific applications have been investigated: (1) The 
measurement of coronary blood flow with the existing CW, small vessel 
flow cuffs, (2) The measurement of coronary blood flow with existing 
pulsed dopplt equipment, and (3) The measurement of aortic flows 
with a catheter- tip, doppler transducer. 

Initial attempts to use the existing CW small- vessel cuffs in the 
measurern • t of coronary velocities have been unsuccessful. A new 
version hits been designed which would be optimized for the measurement 
of coronary flow. This new transducer assembly will be constructed 
and tested in the near future. 

Acute studies have deomons trated the difficulties involved in using 
pulsed flowmet. rs and transducers to measure flow in very small vessels 
such as the coronary arteries. These problems arise partly from the 
relatively long ivelength of the six megahertz ultrasound in compari- 
son to the size of the coronary lumen. Experiments are continuing 
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to determine those combinations of transducers , construction and 
excitation frequency that would facilitate the measurement of coronary 
blood flow by pulsed Doppler flowmeters. 

A concurrent effort has been initiated to combine an SEL directional 
flowmeter and a commercially available, catheter tip, transducer 
produced by Southwest Rescar h Institute, San Anton .'o, Texas. This 
catheter, positioned by use f X-ray techniques, could be used for 
(1) relatively non- invasive measurements of velocities in the major 
blood vessels prior to (or in place of) .implantation of a more com- 
plicated telemetry system, (2) comparison with previously implanted 
electromagnetic flowmeters to establish a non-occlusive "flov. zero” 

(3) comparison with, previously implanted non-directional Doppler flow 
meters to add directional information, and (4) compari on with waveforms 
produced by transcutaneous Doppler flowmeters. 

The first waveforms collected are shown in the figure, below. 
Insinuation vras recorded while the catheter was positioned along the 
axis of flow In the descending aorta. Research on thi f technique 
will continue. 



Velocity in the Descending Aorta — Catheter- tip Doppler 


i 
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III. CONCLUSIONS 

Sample implantable flowmeter designs, suitable for use with a dog' 
aor'a, human aorta, human renal artery and human coronary artery, are 
tabulated at the er of Section A, on the Implantable Puls- d Flowmeter, 
together with a figure of merit for each, indicative of the degree of 
precision in measurement technique possible, and I he way in which this 
is dependent on physiological cons traints . Future work in this area 
will be directed toward confirmation of these designs in actual construc- 
tion and tests. 

For transcutaneous observation of blood flow, several systems have 
been designed. One design concep ' , the H- array, has been constructed 
and tested satisfactorily using target specimens in a water tank. 

Lumen area and vessel oriental ion can be determined by this means. 

This information, combined with Doppler information from the returned 
echoes enables calculation of flow volume. Two simplified new concepts 
ha' e been developed recently, which avoid the need to determine vessel 
orientation. The simpler of these will be constructed during the ne't 
year. It will have a single linear array for visuali ition and a 
separate transducer to produce uniformly distributed sound for use 
in measuring the first moment of the Doppler sj ectrum. Wh-n. the velocity 
distribution is not symmetrical across the lumen, the first moment of 
the Doppler spectrum is the only type of average measurement that 
reflects accurately the tr ; velocity. 

Technique development for implementing flowmeter use resulted in 
success with CW flowmeters but some problems with pulsed flowmeters on 
small vessels. Higher frequency sou.d, than the current 6 MHz, will 
a3 o be needed for satisfactory operation of the pulsed Doppler flow- 
meters on small vessels. Its shorter wavelength should permit more 
accurate observations to be made, but several design criteria must 
be modified as a result. 



